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A search for the decay of the r lepton to five charged and two neutral pions is performed using 
data collected by the BABAR detector at the PEP-II asymmetric-energy e^e~ collider. The analysis 
uses 232 fb~ of data at center-of-mass energies on or near the T(45) resonance. We observe 10 
events with an expected background of 6.51^4 events. In the absence of a signal, we set the limit 



on the branching ratio B(t~ — » ?,-k~ 2-k'^ 2-k" v-r) < 3.4 x 10~ at the 90% confidence level. This 
is a significant improvement over the previously established limit. In addition, we search for the 
decay mode t~ — + 2uj-k~ 
calculate the upper limit B{t~ 
first upper limit for this mode 



We observe 1 event with an expected background of OA'^qI events and 
2LLin~i'r) < 5.4 X 10~^ at the 90% confidence level. This is the 



PACS numbers: 13.35.Dx, 13.66.De, 13.85.Rm 



Hadronic decays of t leptons provide an excellent lab- 
oratory for the study of the strong interaction. Decays 
of the T with one or three charged particles in the fi- 
nal state have been well studied in the past []|. Higher 
multiplicity decays, however, have considerably lower 
branching ratios ll|, and high luminosity experiments 
are needed to study their dynamics and search for new 
modes. The BABAR experiment has recorded a large sam- 
ple of e+e^ -^ T^T~ events suitable for detailed searches 
for high multiplicity r decays. 

The T" — > 37r~27r+27r°i'T- mode J2| is of particular in- 
terest. It has not been observed yet, and only an upper 
limit B{t- -^ 'i'K-2'K+2TT^Vr) < 1.1 X 10"'' at the 90% 
confidence level (CL) has been set by the CLEO collabo- 
ration [3J . The reason for the suppression of seven-pion t 
decays is the limited phase space of this decay [j, y] . For 
the same reason, if this decay is observed with sufficient 
statistics, it may lead to a more stringent limit on the r 
neutrino mass. 

Since r decays to five charged pions and a 7r° meson 
most likely involve resonances (e.g., oj or rj) [Q], it is ex- 
pected that the r" -^ 37r"27r+27r*'i^7. decay may also 
proceed through resonances. According to calculations 
based on isospin symmetry [2|, the decay t~ -^ 2ujT:~v.r 
is expected to be the dominant mode. 

This analysis is based on data recorded with the BABAR 
detector at the PEP-II asymmetric-energy e+e" storage 
ring operated at the Stanford Linear Accelerator Center. 
The data sample consists of 232 fb~^ recorded at center- 
of-mass (CM) energies of 10.58 GeV and 10.54 GeV. With 
an expected cross section for t pairs of (Trr = (0.89 ± 
0.02) nb [S|, the number of produced r-pair events is 
Nrr = (206.5 ±4.7) X 10^ 

The BABAR detector is described in detail in Ref. [9j, 
and only a brief description is given here. Charged- 
particle momenta are measured with a 5-layer double- 
sided silicon vertex tracker (SVT) and a 40-layer drift 
chamber (DCH) inside a solenoidal magnet with a 1.5 T 
magnetic field. A calorimeter (EMC) consisting of 6580 
CsI(Tl) crystals is used to measure the energy of elec- 
trons, positrons, and photons. A ring-imaging Cherenkov 
detector is used to identify charged hadrons, in combina- 



tion with ionization energy loss measurements in the SVT 
and the DCH. Muons are identified by an instrumented 
magnetic- fiux return (IFR). 

Monte Carlo (MC) simulations are used to estimate 
the T~ -^ 37r"27r+27r°z/T- signal efficiency and background 
contamination from other r decay modes. The produc- 
tion of r pairs is simulated with the KK generator jl^ , and 
non-signal t lepton decays are modeled with TAUOLA [llj 
according to measured rates JlJ]. The background pro- 
cesses e'*'e~ -^ qq {q = u,d,s,c,b) are simulated us- 
ing the JetSet package Jl^. Signal events are gener- 
ated using phase space with a, V — A interaction. We 
find no significant variation in efficiency within the phase 
space. The simulation of the BABAR detector is based on 
GEANT4 |T^ . 

The principal backgrounds to our signal come from 
e+e~ —>■ qq processes and multi-pion r decay modes 
involving at least one tt", namely t~ -^ Sn^ 2Tr^ ir'^ Vr , 
T" -^ 27r~7r+27r°i/T- and t~ -^ 27r"7r+37r"i/T- modes. The 
r" — > 37r~27r"'"7r''z/T- contribution comes from reconstruct- 
ing an additional (fake) tt" , while the three-prong modes 
contribute through the tt^ decay to a photon pair and 
subsequent photon conversions in detector material. 

The event selection criteria were developed to suppress 
the background while maintaining high signal efficiency. 
Events with six charged particle tracks and a net charge 
of zero are first selected. To ensure well-reconstructed 
tracks, each track is required to have a minimum trans- 
verse momentum of 100 MeV/c, a distance of closest ap- 
proach to the interaction point in the plane transverse 
to the beam axis (DOCAxy) less than 1.5 cm, and a 
distance of closest approach along the beam direction 
less than 10 cm. Four or more tracks are required to 
have hits in at least 12 DCH layers. Photons are recon- 
structed from clusters in the EMC and are required to 
have a minimum energy of 50 MeV, energy deposited in 
at least three crystals, and a lateral energy profile con- 
sistent with that of a photon. In addition, to suppress 
background from backscattering in the EMC, the angle 
between the position of a cluster and the impact point 
of the nearest charged track at the EMC surface, as seen 
from the interaction point, is required to be more than 



0.08 radians. 

The TT*^ mesons are reconstructed from two photon can- 
didates passing the photon selection criteria described 
above. We first search for tt^ candidates with energy 
E^o > 450 MeV and mass 113 < M^^ < 155 MeV/c^. If 
two or more tt*^ candidates share a photon, only the one 
with the smallest \M^j - M^o^^l, where M^P'^ value is 
taken from (l|, is retained. Next, we repeat the procedure 
for 7r° candidates with energy 300 < Ej^o < 450 MeV and 
mass 120 < M^-y < 148 MeV/c^. 

The T pair is produced approximately back-to-back in 
the e+e^ CM frame. This allows the event to be di- 
vided into two hemispheres by a plane perpendicular to 
the thrust axis, where the thrust is calculated from all 
charged particles and photons in the event [l^. The 
event thrust magnitude is required to be larger than 0.9. 
This requirement rejects more than 90 % of the qq back- 
ground and the e+e^ -^ BB background is suppressed to 
a negligible level. Events are required to have one track 
in one hemisphere (the tag side) and five tracks in the 
other hemisphere (the signal side). To further suppress 



the background from e^ 



qq events, we demand a 



well-identified electron or muon on the tag side with at 
most one additional photon with energy E^ < 500 MeV. 
The combined mass of all charged particles and photons 
in each hemisphere is required to be less than 3 GeV/c^. 
Finally, only events with exactly two 7r° candidates on 
the signal side are kept for further study. The efficiency 
of the two 7r° selection in the signal MC is 13.0 %. 

The visible energy, defined as the sum of the CM 
energy of the charged tracks and the reconstructed 7r° 
mesons, is required to be less than the CM beam energy 
£'boam = 5.29 GeV in each hemisphere of the event. The 
residual energy i?rcs, defined as the neutral energy on the 
signal side not associated with the reconstructed r decay 
products, is required to be less than 300 MeV, reducing 
the background from e^e^ — > qq and r^ -^ 2tt^ tt'^ Stt^ I'r 
events. 

To reconstruct the signal event, an approximation of 
the T invariant mass is used: 
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where the r neutrino is assumed to be massless and travel 
along the direction of the combined momentum vector 
Pjii of the seven hadrons and its energy is taken to be 
the difference between i^beam and the combined energy 
Ejh of the hadrons in the CM system. The variable M* 
is called the t pseudo mass [lj|, and its distributions for 
signal and background MC events are shown in Figure ^ 
The advantage of M* over the invariant mass M^h is 
a considerably better separation of the signal from the 
hadronic qq background. 

We apply particle identification on the signal side, de- 
manding four out of five tracks to be identified as pions 
with high probability, and apply looser identification cri- 
teria to the fifth track. This requirement significantly 
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FIG. 1: Pseudo-mass distribution below 2 GeV/c'^ for signal 
and background MC samples. Signal is plotted assuming 
3{t~ — > Z'k~2-k'^2-k^Vt) = 10"*. For illustrative purposes 
the 1-prong tagging requirements are not imposed here. 



reduces the background from t events with photon con- 
versions and e+e~ — > qq events containing kaons. 

We further suppress photon conversions by requiring 
the invariant mass of each pair of oppositely charged 
tracks to be larger than 5 MeV/c^. In addition, we apply 
cuts on the sums of the two lowest transverse momenta 
and two largest DOCAxy of the tracks on the signal 
> 
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The final event count is performed in the signal region 
1.3 < M* < 1.8 GeV/c^. According to MC studies, the 
signal efficiency after all cuts is (0.66±0.05) %. The error 
is a combination of systematic and statistical uncertain- 
ties. The systematic uncertainty on the signal efficiency 
includes contributions from the reconstruction of charged 
tracks and photons (4.3%), the reconstruction of two 7r° 
mesons (6.6%), and the uncertainty associated with the 
particle identification on the signal and tag sides (1.7%). 
A statistical uncertainty (1.8%) due to limited MC sam- 
ples is added in quadrature to the systematic uncertainty. 

The simulation of r-pair events yields a reliable es- 
timate of their expected background contribution, ver- 
ified by modifying the event selection criteria to sup- 
press the qq background and allow for more r events. 
The largest background is predicted to come from t~ — > 
'i'K~2'K'^T:^Vr decays. For a detailed study, we use an 
MC sample of t~ -^ in^ 2tt^ tt^ Vr events correspond- 
ing to 1900 fb~ of data. The pseudo-mass spectrum of 
the events passing the selection criteria is fitted with a 
"Crystal Ball" probability density function (PDF) [l^. 
In order to determine the shape parameters of this PDF, 
we first fit a larger sample selected without tagging of the 
one-prong side. Using this fixed shape, we then estimate 
the number of t~ -^ 'iTT~2'K~^'K^Vr events within our sig- 
nal region (1.3 < M* < 1.8 GeV/c^) from the MC sam- 



pie with the one-prong tag appHed. We obtain 3.6±0.6 
events, scaled to the luminosity of 232 fb~ , where the 
uncertainty is statistical only (see Figure|21 left). Simply 
counting the number of events in the signal region yields 
3.2 (scaled) MC events. 

The uncertainty of the t^ -^ 3tt^2tt^tt'^i't- background 
estimate is based on the uncertainties of the fitted PDF 
shape parameters, namely, the central value and the 
width, and the correlation between them. The values 
of the PDF shape parameters are randomly generated 
according to their uncertainties expressed in the covari- 
ance matrix, and the resulting PDF is then used to esti- 
mate the number of background events in the signal re- 
gion. The total uncertainty from the fitting (0.6 events, 
16.7%) is added in quadrature with systematic uncer- 
tainties in the reconstruction of the tracks and neutrals, 
particle identification, luminosity and r-pair cross section 
(8.4 %) and the uncertainty in the branching ratio of the 
T- -^ 37r-27r+7r°j/^ decay mode (14.9%). 

An additional background contribution is expected 
from the t^ -^ 27r^7r+27r"i^^ mode. Using an MC sample 
corresponding to 675 fb~ of data we estimate 0.7±0.5 
background events in the signal region from this source. 
The uncertainty is dominated by the MC statistics. Con- 
tributions from other generic r decays are negligible. 
Combining both sources of the r background, we expect 
a total of 4.3 ±1.0 background events in the data. 

For this analysis, a comparison of MC simulation and 
data has shown that the e+e~ — > qq background con- 
tributions cannot reliably be extracted from simulation 
due to difficulties in modeling the fragmentation pro- 
cesses. The shape of the simulated pseudo-mass dis- 
tribution appears to agree with the shape in the data, 
but the overall normalization does not. Therefore, the 
qq background is estimated directly from the data, by 
fitting the data pseudo-mass spectrum with the sum of 
two Gaussians. This PDF is motivated by MC studies, 
which show that the e+e" -^ (uu, dd, ss) and e+e~ -^ cc 
backgrounds have Gaussian pseudo-mass shapes with dif- 
ferent parameters. The double-Gaussian fit to the MC 
pseudo-mass distribution of qq background is shown in 
Figure El (right). 

To extract the qq background in the signal region, we 
subtract the expected r background contribution from 
the data pseudo-mass distribution, and fit the resulting 
histogram in the range 1.8 < M* < 3.3 GeV/c^ with 
a double-Gaussian PDF whose means and sigmas are al- 
lowed to float. To avoid experimenter bias, this fit is per- 
formed "blind" , with the data in the signal region hidden. 
The fit function is then extrapolated below 1.8 GeV/c^ 
and its integral between 1.3 and 1.8 GeV/c^ yields the qq 
background estimate in the data, 2.2 events. 

To calculate the statistical uncertainty of the qq back- 
ground estimate we vary the number of events in each bin 
of the data qq pseudo-mass spectrum above 1.8 GeV/c^ 
according to its Poisson error and refit the resulting his- 
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FIG. 2: Monte Carlo simulated pseudo-mass distributions of 
the T~ -^ 3tv~ 2-jT^ ty" Ut background with a "Crystal Ball" 
shape PDF superimposed (left) and e^e^ -^ qq background 
fitted with the sum of two Gaussians (right). The distribu- 
tions are not normalized to the data luminosity. 



togram for a new estimate. The statistical uncertainty of 
1|^'q events is extracted from the variance of the distri- 
bution of the generated qq background estimates. Varia- 
tions in the functional form of the fit PDF are taken into 
account as a systematic uncertainty of Iq'q events. The 
total uncertainty is calculated by adding the statistical 
and systematic uncertainties in quadrature. Thus, the qq 
background estimate is 2.2lJ;g events. 

To validate the e+e~ —>■ qq background estimate 
method, we apply it to a r-event-free data sample, ob- 
tained by requiring at least 3 photons with energies 
greater than 300 MeV on the tag side not associated with 
a TT*^. This requirement effectively suppresses r events to 
a negligible level and provides a clean qq sample in the 
data. Comparison between the expected and observed qq 
background levels for this sample shows good agreement, 
11.8 predicted background events vs. 12 observed. 

Another cross-check we perform is the branching ra- 
tio measurement of the t~ — > 3tt~ 2tt^ tt'^ lyr decay mode 
using the same selection criteria (except for demanding 
only one 7r° on the signal side instead of two) as described 
above. The measured branching ratio is consistent with 
the Particle Data Group's value 0|. 

Combining the background estimates from r and qq 
events, we calculate a total of 6.5^^4 background events. 
Figure illustrates the final pseudo-mass spectrum of 
the data, along with the expected background PDF. We 
observe 10 events in the signal region and conclude that 
there is no evidence for the r^ -^ 37r^27r+27r'^i^7. decay. 

The upper limit for the r^ -^ 37r^27r+27r'^i/7- decay 
branching ratio is calculated from 



B{t- -^ 37r"27r+27r°zy^) < 
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where XNi,ig„^i is the upper limit on the number of signal 
events at the 90% CL. This number is obtained using 
a limit calculator program [ig that follows the Cousins 
and Highland approach J17l | of incorporating systematic 
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FIG. 3: Pseudo-mass distribution of the data events passing 
the T~ — > 37r~27r^27r Ur selection criteria. The solid curve 
represents the total expected background PDF. The dashed 
curve illustrates the r background contribution. 



uncertainties into the upper limit, using the numbers of 
expected background and observed events, as well as the 
uncertainties on the background, signal efficiency and the 
number of r pairs. We find Xn^^ ^^^ = 9.2 events and 
B{t~ -^ 37r-27r+27rOi/^) < 3.4 xT6~^ at the 90% CL. 
Table in summarizes the results of this analysis. 



TABLE I: Signal efficiency, expected background, observed 
data events, and the upper limit of the t~ -^ 3'K~2n 
decay at the 90 % CL. 



2n°v^ 



Nrr 



(206.5±4.7)xl0^ 



T~ — > 37r~27r+27r°i/T- efficiency 
Expected t'^t~ background 
Expected qq background 
Expected total background 
Observed events 



(0.66±0.05)% 
4.3±1.0 events 

2.2t^-' 
6.5^ 



-1.0 

h2.0 
-1.4 



events 



events 



10 



Bi 



37r-27r+27r"iy^) 



< 3.4 X 10- 



In addition to this inclusive result, we also search for 
the resonant decay mode t~ —>■ 2L0'K~Vr with the sub- 
sequent decay to — ^ tt^tt+tt'^, which is predicted to be 
the main channel for the r~ -^ 37r^27r+27r*'i'T- decay |7|. 
The T~ — > 2un:~i'r mode has a much narrower allowed 
pseudo-mass range (1.7 < M* < 1.8 GeV/c^) due to its 
kinematics. For the same reason, the background level 
is expected to be much smaller. The event selection is 
re-optimized for this analysis. Photons are required to 
have a minimum energy of 50 MeV, energy deposited in 
at least two crystals and a lateral energy profile consis- 
tent with that of a photon. Reconstructed tt" candidates 
must have energies above 200 MeV. The uj resonance is 



reconstructed as a tt+tt" 
ant mass of 0.76 < AIj^+. 

Reconstruction of both lo mesons suppresses the back- 
ground and therefore further selection cuts can be sub- 



combination with an invari- 
< 0.80 GeV/c2. 



stantially loosened to increase the signal efficiency. The 
conversion veto and the i?res cuts are not used. In ad- 
dition, we allow one charged particle of any type on the 
tag side, and only loose pion identification is required on 
the signal side. As a result, the t~ — > 2ujTT~i'r efficiency 
for this selection is (1.53±0.13) %. The uncertainty is a 
combination of systematic and statistical uncertainties, 
as described above for the inclusive r~ -^ 3tt~2tt^2tt'^i/t 
analysis. 

The background is estimated from MC simulation (see 
Figure 0J. As in the inclusive analysis, while there is 
a discrepancy between the data and MC qq yields, the 
shape of the MC qq pseudo-mass spectrum agrees with 
the data. As a result of the study we expect negligible 
qq contribution in the signal region. The uncertainty on 
the qq background estimate is calculated using the same 
technique described for the inclusive t~ -^ Zti^2'k^2-k^Vt- 
analysis. The total expected qq background is O.O^o^o 
events. An additional contribution comes from the t^ -^ 
lo2t:^ TT^ Vr mode. Out of 530 fb^ of MC simulated 
T^ — > lu2it^'k^i/t events, only 1 event is found in the 
signal region. Thus, we expect 0.4J^q 4 events in 232 fb~ 
of data. The uncertainty in the r background estimate 
is calculated as a Poisson error of 1 event at 68 % CL. 
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FIG. 4: Pseudo-mass distributions of the data (points) and 
MC (shaded histograms) events passing the r~ — > 2ijJ-k~Vt 
selection criteria. The dark shaded histogram corresponds to 
the r background, whose level is determined from the simu- 
lation. The light histogram shows the total background, with 
the level of the qq contribution scaled to agree with the data. 
The data signal region below 1.8 GeV/c^ was blinded during 
the background estimation. 



We find 1 event passing the selection criteria in 
232 fb""'^ of data, which is consistent with the expected 
background. We calculate the upper limit of the t~ -^ 
2u;TT~i'r decay branching ratio using the limit calcula- 
tor [lil, which yields Xn,,^„^i = 3.4 events at the 90 % CL. 
The upper limit for the decay, B{t~ -^ 37r~27r+27r''i/T-) < 
5.4 X 10^^, is significantly lower than for the inclusive de- 
cay r^ -^ 3TT~27T~^2Tr'^i'r- TablellTlsummarizes the results 
of the r^ -^ 2uj'K~Vr search. 
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TABLE II: Signal efficiency, expected background, observed 
data events, and the upper limit of the r^ -^ luj-n^v^ decay 
at the 90 % CL. 



N. 



(206.5±4.7)xl0^ 



2(jj7r 



efficiency 



(1.53±0.13)% 



Expected t^t~ background 
Expected qq background 
Expected total background 
Observed events 



-1.0 
-0.4 

1.0 



0.4ta4 events 
events 



0.4Io;4 events 
1 



B{t- 



2a; TT 



< 5.4 X 10" 



In conclusion, we present results of a search for the 
T^ -^ 37r^27r+27r''i^T and r^ — > 2ijjT:^VT- decay modes 
using 232 fb~ of data collected by the BABAR detector. 
No evidence for these decays is found. We calculate 
B{t- -^ 2,TT-2TT+2Ti^Vr) < 3.4 X IQ-*^ at the 90% CL, 
improving the existing experimental limit for this mode 
by more than a factor of 30. The upper limit for the 
decay, B{t'' — > 2w7r~z^T-) < 5.4 x 10^^, is reported here 
for the first time. 
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